Josephson dynamics of a spin-orbit coupled Bose-Einstein condensate 

in a double well potential 
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We investigate the quantum dynamics of an experimentally realized spin-orbit coupled Bose- 
Einstein condensate in a double well potential. The spin-orbit coupling can significantly enhance 
the atomic inter-well tunneling. We find the coexistence of internal and external Josephson effects 
in the system, which are moreover inherently coupled in a complicated form even in the absence of 
interatomic interactions. Moreover, we show that the spin-dependent tunneling between two wells 
can induce a net atomic spin current referred as spin Josephson effects. Such novel spin Josephson 
effects can be observable for realistically experimental conditions. 
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I. INTRODUCTION 

Based on the Berry phase effect [H, [|| and its non- 
Abelian generalization [3| , the creation of synthetic gauge 
fields in neutral atoms by controlling atom-light inter- 
action has attracted great interest in recent theoretical 
studies [34l8j. and has been realized in spinor Bose- 
Einstein condensates (BECs) in the pioneering experi- 
ments of the NIST group , [2(| and also in several sub- 
sequent experiments of other groups [2l| - |23j . The neu- 
tral atoms in the generated effective Abelian and non- 
Abelian gauge fields behave like electrons in an electro- 
magnetic field [22j or electrons with spin-orbit 
(SO) coupling 20]. Different from electrons that are 
fermions, the atoms with the synthetic SO coupling can 
be bosons and typically BECs. This bosonic counterpart 
of the SO coupled materials has no direct analog in solid- 
state systems and thus has received increasing attention 
[24T - I44T ] for different types of SO coupling, different inter- 
nal atomic structures (pseudospin-1/2, spin-1 and spin-2 
bosons, etc.), and different external conditions (homoge- 
nous, trapped and rotated). These theoretical investiga- 
tions focus mainly on the static properties of SO coupled 
BECs and have reveal rich phase diagrams of the ground- 
"25T-I281 [30l Hl| and exotic vortex structures (33r - 
41] . However, to our knowledge, their dynamics 
has been less studied [24], |42H44| , where the SO coupled 
BECs are demonstrated to exhibit interesting relativistic 
dynamics, such as analogs of self-localization [Hj], Zit- 
terbewegung [43[ and Klein tunneling [44| under certain 
conditions. 

On the other hand, quantum dynamics of a BEC in 
a double well potential has been widely investigated. In 
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FIG. 1: (Color online) A schematic representation of (a) a SO 
coupled BEC in a double well trap, (b) The dynamic process 
of the system, where the blue solid and dashed lines repre- 
sent the inter-well tunneling without and with spin-flipping, 
respectively, the green circles represent the Raman coupling. 
The inter-well spin-flipping tunneling induced by the Raman 
coupling is negligible under current experiment conditions 
[20I . IHol ] . The atomic interaction terms are not shown in the 
figure. 



particular, the coherent atomic tunneling between two 
wells results in oscillatory exchange of the BEC, which 
is analogous to the Josephson effects (JEs) for neutral 
atoms [4oT - |48| . The weakly interacting BECs provide a 
further context [46l - |48| for JEs in superconductor sys- 
tems because they display a nonlinear generalization of 
typical d.c. and a.c. JEs and macroscopic quantum self- 
trapping, all of which have been observed in experiments 
49-511. Apart from the conventional single-species BECs 
45 - 48 j , the Josephson dynamics of two-species BECs [HH 
and spinor BECs without SO coupling [53[ have also been 
studied [5~4T - l57j ] ; however, the dynamics of SO coupled 
BECs is yet to be explored. 

In this paper, we investigate the dynamics of a specific 
SO coupled BEC, which was realized in the experiment 
of the NIST group, in a double well trapping potential. 
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We find that the SO coupling in the system contributes 
to and increases the atomic tunneling to a large extent, 
which can significantly enhance the atomic JEs. The full 
dynamics of the system contains both internal and exter- 
nal JEs, which are moreover inherently coupled in a com- 
plicated form even in the absence of interatomic interac- 
tions. We further demonstrate that the spin-dependent 
Josephson tunneling can lead to a net atomic spin current 
by varying conditions, which we refer to as spin Joseph- 
son effects. The predicted spin- Josephson currents are 
robust against the parameter adjustment and varying ini- 
tial conditions, and can be observable in the SO coupled 
BECs under realistic experimental conditions. 

The paper is organized as follows. In Sec. II we con- 
struct a model that can be used to study the quantum 
dynamics of a SO coupled BEC in a double well poten- 
tial. Then, in Sec. Ill, the Josephson dynamics of the 
constructed system is investigated, with the complicat- 
edly coupled internal and external JEs being addressed. 
In Sec. IV we demonstrate that the spin JEs exhibit in 
the system under realistic conditions. Brief discussions 
and a short conclusion are given in Sec. V. 



II. MODEL 

In a very recent experiment, the NIST group realized 
a synthetic SO coupling in the 87 Rb BEC, in which a 
pair of Raman lasers generate a momentum-sensitive cou- 
pling between two internal atomic states [20| . In the bare 
pseudo-spin basis | f)& = \ m F = 0) and | 4-)& = \nip = 
— 1), the SO coupling is described by the single particle 
Hamiltonian given by [20| 
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(1) 



where p is the atomic momentum in the xy plane, m is 
the atomic mass, S is the tunable detuning behaved as a 
Zeeman filed, k^ is the wave number of the Raman laser, 
and ft is the Raman coupling strength. Such kind of 
SO coupling is equivalent to that of an electronic system 
with equal contribution from Rashba and Dresselhaus SO 
coupling, and thus it is effective just in one-dimension 
(ID). So we restrict our discussions in ID and focus on 
the motion of atoms along x axis by freezing their y and 
z degrees of freedom. 

To proceed further, we introduce the dressed pseudo- 
spins I |) = e~ ikLX \ t)& and | |) = e lkLX \ 4-) b @,|il|, then 
the single-particle Hamiltonian in ID (along x axis) can 
be written as 

ha = — — - + 2ak x a z + -a x + -a z , (2) 
2m 2 2 

where k x is the atomic wave vector operator, and a — 
E r /kL is the SO coupling strength with E r = Ti 2 k\/2m 
being the single-photon recoil energy. The dispersion re- 
lation of the single particle Hamiltonian @ with 6 = 



is E ± (k x ) = ^ ± A/4a 2 fc2 + ft2/4, which exhibits a 
structure of two branches. We are interested in the 
lower energy one E^{k x ). There is only one minimum 
in k x = for large Raman coupling ft > 4E r , where 
the atoms of both atomic levels condense. However, 
the lower branch for ft < AE r presents two minima for 
condensation of dressed pseudo-spin-up (left one) atoms 
and dressed pseudo-spin-down (right one) atoms, respec- 
tively. The Raman coupling and a small 5 modulate the 
population of atoms in these two states [2(| . Here we fo- 
cus on the later regime, i.e. ft < 4E r , because such BEC 
with spin-separated and non-zero central momentum is 
more interesting in contrast to a regular BEC with zero 
central momentum. 

To be clearer, we can rewrite the Hamiltonian @ as 



i _ ( H t ft/2 \ 
h ° { ft/2 H X ) 



(3) 



where H t = £-(k 2 x + 2k L k x ) + f and H x = £-(k 2 x 
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Since it is more straightforward to describe 



the system in terms of dressed pseudo-spin states com- 
pared with using bare ones, we will work in the dressed 
pseudo-spin space and simply refer to dressed pseudo- 
spin as spin for convenience hereafter. We also note that 
the parameters k^ , ft and S in the single-particle Hamil- 
tonian can be tuned independently in a wide range [20| . 
making the SO coupled BEC a suitable platform for in- 
vestigating the Josephson dynamics in the presence of SO 
coupling. 

Now we turn to consider such a SO coupled BEC in a 
double well potential denoted by V(x) as shown in Fig. 
1(a). Note that the double well potential here is assumed 
to be spin-independent. To investigate the dynamics of 
the system, we adopt the two-mode approximation (45I — 
HH with the field operator 



(4) 



where ipj a (x) is the ground state wave function of the j 
well (j = L,R) with spin a (a =j, J,), and a,j„ is the an- 
nihilation operator for spin a in the j well, satisfying the 



bosonic commutation relationship [a 



The validity of the two-mode approximation holds under 
two conditions: the weak atomic interaction and small ef- 
fective Zeeman splitting, as the atoms can not be pumped 
out of the lowest state of each well in this case. In the sec- 
ond quantization formalism, the total Hamiltonian reads 



n 



dx&(x) h + V(x) + h int 



(5) 



where the two-component field operator 'f = (^,^) T , 
and the interaction Hamiltonian hi n t will be specified 
below. Substituting Eq. (0} into Eq. ([5]) , one can rewrite 
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the total Hamiltonian as 



H = gjo-a jo-Qj-ff + (j<ro->aL a aRa< + h.c 
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(6) 



where £j CT = / dxif>j a (x) [H a + V(x)]i/jj a (x) « jfiwj — 
EV is the single-particle ground state energy in the 
j well with uij being the harmonic frequency of this 
well, J aa = J dxtp* La (x)[H a + V(x)]iJj Ra -(x) and J aS = 
J dxip La (x)^ip Ra(x) (with a and a referring to different 
spins) are the tunneling terms shown in Fig. 1(b). In 
addition, the interaction Hamiltonian is given by 
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(7) 



2h J a„ 



ml i - J dx\^j a [x)\ 2 \^jj a/ {x)\ 2 is the effec- 
tive ID interacting strength with a aa i being the s-wave 
scattering length between spin a and a' and l± being the 
oscillator length associated to a harmonic vertical con- 
finement. Note that here we have ignored the inter- well 
atomic interactions because the s-wave scattering length 
(which is on the order of nanometers) is much smaller 
than the inter-well distance (which is on the order of 
micrometers). We have also dropped the inter- well cou- 
pling since its strength is exponentially smaller than the 
intra- well counterpart. The Hamiltonian © describes 
the dynamic process of the system schematically shown 
in Fig. 1(b). 

For simplicity, we assume the double well potential to 
be symmetric as shown in Fig. 1(a), with each well hav- 
ing the same harmonic trapping frequency lo. Thus we 



have El = £r and g t 



(L) 



9, 



(Ji) 



Such kind of double well 



potential can be generated in experiments [51| with the 
form 
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V(x) = a(x 2 - b 2 ) 



(8) 



where the parameters a and b are both tunable in the 
experiments [Hi]]. Expanding V(x) near x — ±6, one ob- 
tains its harmonic form as V^(x) |mw 2 (i ± b) 2 and 
thus a = mui 2 /8b 2 . The ground state wavefunctions of 
the BEC in each well potential with each spin can be ap- 
proximately represented by its corresponding lowest en- 
ergy single-particle wavefunction, which can be worked 
out by solving the equations [H a + ^muj 2 (x ± b) 2 ]tpj a = 
ejaifjja (here ± are for j — L,R, respectively). The re- 
sults are [H 



V>lt = (Po L \x)e- lkLX , 
iJLi = ^ L) (x)e lk ^, 
1 p Rt = l f { R) (x)e- ik ^, 



(9) 
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FIG. 2: (Color online) The energy scales of tunneling terms 
J^ 1 ' as a function of b for (a) lo = 1 /Ltm, (b) lo = 2 /im, (c) 
lo — 3 /jm, and (d) lo — 4 /jm, respectively. In (a)-(d) we set 

n/n = E r /n = 22.5 kHz. 



where ip (x) 



and ifQ^ (x) 



. 1 e ( x ^ l 2l ° with Iq — ^/h/muj being the oscillator 

length. Substituting Eq. © into the expressions of J atT i , 
one can obtain 



Ju = J (T) 



J 



n 



j(T) + j(SO) + j(V) + J (Z) j 
J(SO) + J (V)_j(Z)^ 



(10) 



where the terms = 



■Jdxtp^tp^, and J( fl ) 



j(SO) = 

= 



f dx<pi L) V(x)<pi R \ 

lfdx4 L) 4 R) , and J(«) = qjdx<pi L) e- 2ik '-*<pV t) . 

Compared with the atomic tunneling of a regular BEC, 
the SO coupled BEC in this system exhibits two addi- 
tional tunneling channels, the SO coupling induced tun- 
neling term J^ ^ and the Raman coupling induced one 
. To clarify the effects of these terms in the tunneling 
processes, we need to work out and to compare the energy 
scales of all the terms J (v \ where r\ = {T, V, Z, SO, R}. 
Substituting Eqs. © and ([5J into Eq. (ITU1) . we can ob- 
tain the following analytical solutions 



where £ T = T^riik 

^ 2m ; v 2 



b 2 l 2 



(11) 



iz = 6/2, £so = -E r , and 



(-1 

k L l o. Since the Zee- 



man filed S is independently tunable to the double-well 
structure and should be small, we here further assume 
S <C E r and thus we focus on the comparison among 
j(SO,r,t,V) _ e g- ec |- s £ Zeeman-splitting induced tun- 
neling will be specified in the Sec. IV. 

For lo ~ \/2b, we have £t ~ and £y ~ I^F' anc ^ 
thus l j0 = 8tt 2 (^) 2 > 100. Here we have assumed 
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the same wavelength Al = 2-K/ki, = 0.8 /xm and recoil 
frequency E r /h = 22.5 kHz as those in the experiments 
[20I ] , b and Iq to be on the order of micrometers [13, [E3 ■ 
In fact, in the regime of 6 2 /Zq ~ [0.5, 2], we find that 



so 



> 



100max{|fr|,|^|}. 



(12) 



Besides, one can check that £r ~ for Zo ~ 1/xm 

and Ai = 0.8/im, and thus the Raman-coupling induced 
tunneling is negligible in this system. The comparisons 
among for some typical parameters are shown in Fig. 
2. In other words, we find that under realistic experiment 
conditions [2Cj, [5CJ, the spin- flipping tunneling induced 
by Raman coupling is negligible but the SO coupling in- 
duced tunneling term j( so ' dominates and moreover it 
greatly enhances atomic tunneling in this system. Thus 
we may rewrite the tunneling terms as 



Jit « 0, 



(13) 



where 7 = exp(— b 2 ~ [0.1,0.6]. It is worthwhile to 
note that the new tunneling terms J^ so ^ and in 
this SO coupled system are both tunable, enabling us to 
study the interesting effects of SO coupling in the atomic 
inter-well tunneling. For instance, one can decrease the 
effective wave number in x axis to the scale fcz, ~ 1/io so 
that ~ 0.37f2 and then the Raman-coupling induced 
tunneling can revive. This can be achieved by adjusting 
the angle between the applying Raman lasers and the 
trapping potential or alternatively by using lasers with 
larger wavelength. In addition, in the same way one can 
tune the recoil energy E r to identify the enhancement 
of atomic tunneling due to the SO coupling (i.e. the 
effect of j( s °}) in experiments. As a first step to inves- 
tigate the system under current experiment conditions, 
we here concentrate on the tunneling regime governed by 
Eq. (USD. 



Hint — in Eq. (J6)), in which the single-particle Hamil- 
tonian is given by 

n a ~ J tt (a[ t a Rt + a LT a^ t ) + J u (a^a^ + az+a^J 

(14) 

Here we have dropped the tunneling terms J aa since 
these spin-flipping tunneling progresses can be negligible 
in the current experiment conditions (20l . 150] ] . In order 
to study the dynamic properties of the system, we need 
to work with the equation of motion. The corresponding 
Hcisenberg equations read 



ifi—a- 



[aj a ,%o] — Juo-Clja + §Ojct + (— 1) 



2 j c ' 



(15) 

where a and a refer to different spin, while j and j 
to different wells, and p = 0, 1 are for =t,|, re- 
spectively Using the mean-field approximation, one has 
hja — (aja} = cLja with djo being c numbers. Thus we 
can rewrite the equations of motion as 

iM ja = Jaadja + § + (~l) f % a ]o- (16) 



By defining a four-component wavefunction $ = 
{a.L^,aLi,aB,t,aRi) T , Eq. (JTSJ is rewritten as ifi^ = 
Hm&, where the Hamiltonian of the system is given by 



Hm = 



2 
2 



n 7 

2 J n 



o I 



\ 

n 

2 . 



(17) 



We now look into the JEs in this system. Let us fur- 



ther express a ?(T as a ?(T 



r , where the particle 



numbers Nj a and phases 9j a are all time-dependent in 
general. According to Eq. ([To]) , we can obtain 



III. FULL DYNAMICS OF THE SYSTEM 

We are now in the position to investigate the quan- 
tum dynamics of the system constructed in the previous 
section. We first address the non-interacting case, i.e. 



ih^f- - hNj J ja = J aa y/Nj a N- ja ^'- e ^ (lg) 
+%^N~N^e i ^-e i .) + (-l)ffjV^. 

Separating the image and real parts of Eq. (fT8|) yields 
two groups of equations as 




n /Ni S 



• cos 



'jo- 



1 ) 2ft- 



(19) 



Eq. (|19p consists actually of eight coupled equations. To 
simplify these equations, we introduce (f> a — 0R a — 9 La 
and p a = Ar ct — for the phase and particle number 
differences between two wells with the same spin er, and 



4>j = — 9j^ and pj = Nj^ — for the phase and 
particle number differences between two spins in the same 
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well j, respectively. Thus we can obtain 

p a = Ci singer + - y^^-l) q £ 2 sin <f>j, 
j 

Pi = ^(-l) p Asm</v + £ 2 sm^, (20) 



where L x = -^^N Ra N La , C 2 = -f-^/N^N^, and 
q = 0, 1 for j = respectively (p = 0, 1 for er =t,|, 

respectively). From the above Eq. (|20|) . we find the 
coexistence of internal JE related to Pj{4>j) and external 
JE related to p Moreover, the internal and external 

JEs are inherently coupled in a more complicated form. 

Before ending this section, we briefly discuss the 
weakly interacting cases, which have been assumed to 
meet the requirement of two-mode approximation. In 
this regime, the mean-field analysis still works well, and 
the dropped term Hint can be taken into count within 
the previous discussions. This leads to two additional 
terms related to interactions into Eq. (|16|) . and now the 
equations of motion are given by 

iTlQjja — JfjfjClja ~\~ ~^^ja ~1~ ( 1)^ 1 '2® J 3' 7 (21) 
~~\~9<7<7 I t^jcr | (Ljcr ~t" Qoo I^j'ct I Qjcri 

where the interacting strength g aa * can be found as 
q rsn t = ^l 1 a .?°' . The estimation of the interaction en- 



By defining the normalized inter- well particle number dif- 
ference for spin a as Z a = [Nr^ — NL a \/N t (—1 < Z a < 
1), the equations of motion (I19p become rather simple 
in this case (similar to those for the regular two species 
BECs 54]), which are given by 

Za = -^ifVl-^sin <f> a , 

j>° = V 7=? cos^ CT + ^z a + ^z s + 

(23) 

The spin-dependent atomic density current is given by 

I a =NfZ a . (24) 
From Eq. (|24[) . we can define the net spin current as 

h =If~h, (25) 
and the total atomic current as 

Ia = h+h- (26) 

We first consider the JEs of the system in the noninter- 
acting limit, i.e. U aa — U aa — in Eq. ([2"3")l . which can 



ergy and the Josephson dynamics in the presence of weak 
interactions will be presented in the next section. 



IV. JOSEPHSON EFFECTS IN WEAK RAMAN 
COUPLING REGIMES 

In the preceding section, we have shown that the SO 
coupled BEC in a double well potential exhibits the com- 
plicated coupled external and internal Josephson dynam- 
ics. We, in this section, consider a specific dynamic pro- 
cess of the system in the weak Raman coupling regime 
(i.e. Q,/E r <§; 1), where the external Josephson dynamic 
dominates. In fact, the manipulation and detection of the 
SO coupled BECs in this regime have been performed in 
experiments [20j |. 

For the weak Raman coupling, we find that the ra- 
tio v = iJo-o-l/f] can reach several hundreds from Eq. 
(fT3"|). Thus within the time scale r ~ Ti/Q, ~ 45 ms 
for Q, = 0.001i? r , one can ignore the effects of the spin- 
flipping tunneling, which leads to two external Josephson 
tunneling processes for different spins. The spins in this 
regime are conserved and then the total particle number 
of spin a N at = + are time-independent con- 
stants. We assume N-ft — N±t — Nt for simplicity. The 
equations of motion (|21[) in this case can be rewritten as 



(22) 

I 

be realized by Feshbach resonance [EH . Under this condi- 
tion, the two external Josephson tunneling processes for 
different spins are decoupled. We numerically calculate 
Eqs. (|23|) . with some typical results of the time evolu- 
tion of Z a for different initial conditions being shown in 
Fig. 3. In the calculations, we have assumed the zero 
Zeeman filed 5 = in Fig. 3(a) and (b), and small Zee- 
man field 5 = 0.01EV in (c) and (d). Compared with 
zero Zccman field small Zeeman field results in 

a deviation in Josephson tunneling strengths J aa and in 
time-cumulative phases 8/2% for different spins. Here 
Zf(t) and Z±(t) demonstrate the oscillatory Josephson 
tunnelings which are similar to the early results in Ref. 
[47| . As shown in Fig. 3(a-d), they are spin-dependent 
and the dynamic evolution of each one depends on its 
own tunneling strength, phase and initial conditions. We 
also calculate I a , I s and I a in this regime with typical 
results being shown in Fig. 4. It is interesting to see 
that the spin-dependent atomic density currents due to 
the spin-related Josephson tunnelings give rise to a net 
spin current (cf. Fig. 4), and moreover in some certain 
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FIG. 3: (Color online) The time evolution of Z a in nonin- 
teracting limitation. In (a) and (b) we have 5 — and 
Jft ~ J-U = — 0.lE r ; In (c) and (d) we have S — 0.01iJ r , 
J-ft = — 0.905i?r, and J4.4. = —0.105-E r . The initial conditions 
are Z t (0) = -Z±(0) = 0.3, t (O) = 0.5^(0) = tt/4 in (a) 
and (c); and Z t (0) = Z±(0) = 0, c/> t (0) = ^(0) = tt/4 in (b) 
and (d). 



FIG. 5: (Color online) The time evolution of Z a in the weakly 
interacting regime with U aa = 0.01-E r and U aa i = 0.01 1E T . 
Other parameters and initial conditions in (a)-(d) are the 
same with those in Fig. 3(a)-(d), respectively. 




2 4 6 8 10 

t / ms 



FIG. 4: (color online) Josephson currents in the noninter- 
acting limit. The time evolution of spin-dependent atomic 
currents I a , a net spin current I s and a total atomic current 
I a in (a) for the same conditions in Fig. 3(a); and in (b) for 
the some conditions in Fig. 3(b). 



initial conditions the total atomic current can be zero, 
which leads to a new interesting pure spin currents (cf. 
Fig. 4(b)). We call such new JEs as spin Josephson 
effects, which can be observable in experiments by mea- 
suring the time-evolution of spin-dependent population 
imbalance of the atomic gas [59(. 

For weakly interacting cases, we have to estimate 



the interaction energy U arT and U aS , which should be 
Uaa 1 U a9 "C Tilo due to the two- mode approximation. 
This requirement results in U aa s» U aS <C 0.lE r . In 
this regime, the two spin- Josephson tunneling processes 
are coupled via atomic interactions. To understand the 
effects of the interaction, we show in Fig. 5 some typical 
results of the time evolution of Z a for the same initial 
conditions and parameters in Fig. 3. It clearly demon- 
strate that the modification of Z a (t) due to atomic in- 
teractions is not significant and even very minor in some 
cases (such as the case for Z CT (0) = and S = in Fig. 
3(b), 5(b)) since the interaction energy is small compared 
with the tunneling energy. Therefore the spin Josephson 
dynamics still exhibit a similar oscillatory feature in this 
regime. 

To see more clearly the oscillatory properties of the 
spin JEs, we have numerically calculated the frequency 
spectra of the net spin currents I s for various conditions 
(such as those in Fig. 3 and Fig. 5). We find that the 
spectra for different cases exhibit a single peak centered 
at the slightly shifted frequency, as seen in Fig. 6. The 
single-peak feature shown in Fig. 6 implies that the spin 
current I s (t) can well be described by a sin- function, 
while the weak interatomic interactions or the small Zee- 
man field can merely modify the period and amplitude 
of the current slightly. Thus we conclude that the spin 
JEs in this system are robust against the parameter ad- 
justment and initial conditions. 
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FIG. 6: (Color online) Spectra of the net spin currents I 3 for 
the cases in Fig. 3(a,c) and in Fig. 5 (a,c). A single peak in 
the spectrum of each case implies that the spin current I s (t) 
is well described by a sin-function. 



atomic inter- well tunneling and consider only the internal 
dynamics in each single well. The atomic tunneling be- 
tween two spins refers to spin-flipping is induced by the 
Raman coupling, and thus such Josephson tunneling is 
in the spin space. Considering the atomic gas condenses 
with a finite but opposite momentum for different spins, 
the internal JEs connect interesting quantum tunneling 
in the momentum space. 

In summary, we have investigated the quantum dy- 
namics of a SO coupled BEC in a symmetric double 
well potential. The SO coupling contributes to atomic 
tunneling between wells and significantly enhances JEs 
for realistic conditions. We have predicted a novel spin 
Josephson effects which can be observed in a practical ex- 
periment since all the required ingredients, including the 
SO coupled BECs and the tunable double well potential, 
have already been achieved in the previous experiments. 



V. DISCUSSION AND CONCLUSION 

Before concluding, we briefly discuss another specific 
dynamic process of the system in the relatively strong 
Raman coupling regime, 3> |</<t<t|, which can be real- 
ized such as by tuning J aa ~ — 0.lE r and f2 ~ E r . In this 
regime, within the time scale h/\ J acr \ one can ignore the 
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